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Abstract. We exactly diagonalize the finite-size XY model with periodic boundary 
conditions and analytically determine the ground state energy. We show that there are 
two types of fermions: singles and pairs, whose dispersion relations have a completely 
arbitrary gauge-dependent sign. It follows that the ground state is determined by a 
competition between the vacuum states (with a suitable gauge) of two parity sectors. 
We finally exhibit some points in finite size systems that forerun eriticality. They are 
associated to single Bogoliubov fermions and to the level crossings between physical 
and unphysical states. In the thermodynamic limit they approach the ground state 
and build up singularities at logarithmic rates. 



1. Introduction 

The analysis of one dimensional spin chains is a useful approach to the modeling 
of quantum computers [T]. This class of systems has been deeply studied in the 
thermodynamic limit [21 [3l H] ; however, experimental and theoretical difficulties impose 
strong bounds on the realization of large scale systems, and this has boosted a high 
interest in finite size systems O El d E]. The investigation of the last few years 
has focused on entanglement [9], [10] in diverse finite-size models, by means of direct 
diagonalization [TT], [121 CSl [HI [IS, [IE]. These studies were boosted by the recent 
discovery that entanglement can detect the presence of quantum phase transitions 
[ni[l8l[l9l[2ni[2l]. 

In this article we exactly diagonalize the XY model with periodic boundary 
conditions, describing a one dimensional chain made up of a finite number of two level 
systems (|-spins) with nearest neighbors coupling, in a constant and uniform magnetic 
field. The XY model is a class of Hamiltonians distinguished by a different value of the 
anisotropy coefficient, which introduces a different coupling between the x and the y 
components of the spins (in particular the isotropic case, corresponding to the case in 
which the anisotropy coefficient vanishes, is known as XX model). 

As for infinite chains [1], the diagonalization procedure is divided in three steps: 
the Jordan- Wigner transformation, a deformed Fourier transform (generalizing the 
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discrete Fourier transform), and a gauge dependent Bogoliubov transformation. After 
the Jordan- Wigner transformation the Hamiltonian, expressed as a quadratic form of 
annihilation and creation operators of spinless fermions, is characterized by the presence 
of a boundary term [2] whose contribution, which scales like 0{1/N) in the calculation 
of real physical quantities, cannot be neglected for finite size systems. However, this 
boundary term vanishes in Fourier space if the discrete Fourier transform is deformed 
with a local gauge coefficient, depending on the parity of the spins anti-parallel to the 
magnetic field [22] . 

There will emerge two classes of fermions, coupled and single ones (in particular 
for the XX model there are only single fermions). The last step of the diagonalization 
procedure is the unitary Bogoliubov transformation, given by a continuous rotation for 
fermion pairs and by a discrete one for single fermions. We will show that this unitary 
transformation is gauge dependent, since it is given by two possible continuous rotations 
for fermion pairs and by either the identity or the charge conjugation operator for single 
fermions. From this it follows that the sign of the dispersion relation is completely 
arbitrary, apart from the constraint that fermions belonging to the same pair have the 
same sign. 

From the arbitrariness of the Bogoliubov transformation it follows that a possible 
expression for the diagonalized Hamiltonian is such that for successive intervals of the 
magnetic field the vacuum energies of the two parity sectors alternatively coincide 
with the ground state and the first excited level: we will exhibit this mechanism of 
"competition" between vacua. 

Finally we show that in finite size systems one can find the "forerunners" of the 
points of quantum phase transition of the thermodynamic systems. They are associated 
to single Bogoliubov fermions and arise at the level crossings between physical and 
unphysical states. At the values of the magnetic field corresponding to the forerunners 
the second derivative of the ground state energy scales as — log A^. Since in the XX 
model all Bogoliubov fermions are single, one re-obtains the well known result that 
in the thermodynamic limit the anisotropic case presents two discrete quantum phase 
transitions whereas the isotropic or XX model is characterized by a continuous one [1]. 

2. The XY Hamiltonian 

We consider A^ spins on a circle with nearest neighbors interaction in the xy plane 
and with a constant and uniform magnetic field along the 2;-axis. The Hilbert space is 
Ti. = <S)ieZN where i)i = is the Hilbert space of a single spin, and Z^v, labeling 
the positions on the circle, is the ring of integers modA^ with the standard modular 
addition and multiplication. The XY Hamiltonian is given by 

H,i9) = -J Yl 
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with 

a- = 1 (g) 1 (g) . . . (g)cr' (g) . . . (g) 1, 2 e Ztv, le{x,y,z} (2) 
where cr' acts on the i-th spin and may be represented by the Pauh matrices, 

J > is a constant with dimensions of energy and G M and 7 G [0, 1] are two 
dimensionless parameters: the first one is proportional to the transverse magnetic field 
and the second one is the anisotropy coefficient and denotes the degree of anisotropy 
in the xy plane, varying from (XX or isotropic model) to 1 (Ising model). As is 
well known, in the thermodynamic limit, the diagonalization of the XY Hamiltonian 
is achieved by means of three transformations: the Jordan- Wigner (JW), Fourier and 
Bogoliubov (BGV) transformations. We will analyze in detail how the topology of the 
circle will induce a deformation on these transformations in finite size chains. 



2.1. Jordan- Wigner and deformed Fourier transformations 

The Jordan- Wigner transformation is based on the observation that there exists a 
unitary mapping 

between the Hilbert space of a system of spins 7i = (C^)*^^ and the fermion Fock 
space JF_(C^) of spinless fermions on sites. Here, 

-^-(f))=Q-0[)", (5) 

n>0 

where f)" = f)*^" for > 1, f)° = C, and Q- is the projection onto the subspace of 
antisymmetric wave-functions [23]. In order to simplify the notation, in the following 
we will use the above isomorphism and will identify the two spaces Ti. = JF_(C^) without 
making no longer mention to U. By virtue of this identification we can consider the 
canonical annihilation and creation JW fermion operators 



^1=( n = e'^"^%+ VzGZ^, (66) 

where af" = (af ± icrf )/2 and riii is the operator counting the number of holes (or spins 
down) to the left of i 

j&N,j<i 

Note that the above definitions rely upon the following (arbitrary) ordering of Z^r: 
[0] < [1] < ... < [A^ — 1], where [k] = k + NZ. In particular, if the choice of the 
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successive elements can be considered natural, and is well adapted to the Hamiltonian 
([1]), the choice of the first element [0] is totally arbitrary and is related to the choice of 
a privileged point of the circle. 

The JW operators anti-commute both on site and on different sites (see (l8cl) ) 
whereas the Pauli operators anti-commute only on the same site: Vi, j G Zat 

{q,Cj} = 0, (8a) 

{4,c]} = o, m 

{cuc\] = 5ij (8c) 



and 

{cr±a±} = for^ = j, [a± a±] = for z ^ j. (9a) 

From equations ( l6aj) -( l66j) one sees that the terms in the Hamiltonian describing the 
coupling between spins [0] and [A^ — 1] = [—1], when written by means of the JW 
operators, is characterized by an operator phase, at variance with the other coupling 
terms; for example the terms coupling the spins along the x axis become 

= Cjc]+i + CjCj+i + c]+ic] + Cj+ic] , Vj e Zjv\{[-1]}, (10a) 

^hH^fo] = e^"("i+^) (c[_i]cfoj + C[_i]C[o] + cfo]cf_,j + C[o]c|_,]) , (106) 

where the number operator, 

n^= ^(l-c]c,), (11) 

counts the total number of spins down in the chain. This introduces some difficulties in 
the diagonalization of the Hamiltonian, because its expression written in terms of the 
fermion operators is characterized by the presence of a boundary term with the same 
operator phase found in equation ( 1106|) 

H^{9) = -J\Y^ [^(1 - 2cic]) + Cjc]+i + Cj+ic] + 7(ciCj+i + cj+^cj) 

- (e^""^ + 1) [(q-i]4o] + ^ni-i]) + 7(q-i]qo] + clf\_,{)\ |. (12) 

In the thermodynamic limit the boundary term can be neglected since it introduces 
corrections of order 1/A^; the problem is then reduced to the diagonalization of the so 
called "c-cyclic" Hamiltonian [2] and can be easily achieved by means of the discrete 
Fourier transform 

Cfc = e"^c,-, VA; G Zn- (13) 

Since we are interested in finite size systems, with finite iV, the boundary term cannot be 
neglected. The main difficulty introduced by the boundary term in the Hamiltonian (fT2ll 
is that it breaks the periodicity of the JW operators, due to the arbitrary dependence 
of the phase e"^"''^ on the ordering of the spins on the circle. This phase clearly depends 
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on the state the Hamiltonian is apphed to. However, equation (I12p can be simphfied 
by noting that the parity of the number of spins down, 

P = e*^"i, (14) 

is conserved 

[P,if^]=0, (15) 
ahhough not so the spin-down number operator n| itself. Its spectral decomposition is 

V=Y,QPo = P+-P-. (16) 

Q=±l 

where 

p+= 

n| even 
■«! odd 

are the projection operators belonging to the eigenvalues q = ±1 olV respectively, and 
|?7,j) is the eigenstate of rij with eigenvalue nj^. Since parity is conserved (equation f|T5l) ) 
the Hamiltonian can be decomposed as 

= P+H^P+ + P^H^P- = H^+^ + H\-\ (18) 

and the analysis can be separately performed in each parity sector, where V acts as a 
superselection charge. 

In each sector the XY Hamiltonian can be diagonalized by deforming the discrete 
Fourier transform by means of a local gauge aj {j G Zat), 

The inverse formula reads 

" ~7n Y^^^KW ^ ^ ^ ^^^^ 

V-'' ucir \ / 



This deformation preserves the anti-commutation relations in the Fourier space 

{cfc,Cfc'} = 0, (21a) 

{44} = 0, (216) 
{cfc,cj.,} = (21c) 

V/c, /c' G Zjv- The local gauge exp (27riaj/A^) can be determined by imposing that the 
Fourier transforms of fllOap and (11061) have the same form. Considering the first terms 
in the sums one gets 

c,cj+i =e'^("^-"^+^)^ e'^b-'^-^^'+^^^'lcfccl,, (22a) 

e'-("^+i)c[„i]cfoj=e^^("i+^)eT(°[-ii-°[oi)l ^ e^[(-i)^-°'='lcfcci„ (226) 
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where e^'^'^"^"'"^^ is uniquely defined in the sector; they have the same form when, Vj G Zjv, 
exp (^("j - "i+i)^ = exp (i7r(n| + 1)) exp (^^("[-i] - "[o])^ • (23) 

Therefore, the left hand side, like the right hand side, must not depend on j: 

aj+i - aj = a, Mj E Z^, (24) 

with a solution to the equation 

exp (27ria;) = exp (ivr(nj^ + 1)) , (25) 

and the phase associated to the first site ao completely free. The solutions in the two 
parity sectors are 

. r O(modiV) iig=-l {n^ odd) 

a = ^—^ (modN) = { 1 (26) 
4 1^ - (modA^) a Q = +1 (nj even). 

Summarizing, by substituting the (sector dependent) deformed Fourier transform 

Cj = ^y=- ^ exp ( ^Tp(^ + a) ) Cfc, 3 e ^jv, (27) 
into equation f|T2l) we obtain 

H^'\9) = - ^ E + 22^4 [cos (2n^) - 

27r^^ j (^eTcgCfc + q- — c\c\^ \ P„ (28) 

where 'ik eI^n 

k = -2a -k + NZ. (29) 

A comment is now in order. Note that, alternatively, instead of the Fourier transform 
one could have deformed the JW transformation in the following way 

cj = e^^^ie-^^^'^+^o^^, Vj G Zn, (30) 

and would have obtained the same results. 

2.2. The Bogoliuhov transformation 

Observe that when 7 > the last term in Hamiltonian (|28|1 couples fermions with 
momenta k and k. In fact, there are two types of fermions, the single and the coupled 
ones (fermion pairs). Their momenta k belong to the two sets 

S,= ^kEZN k = k'^ = i^kEZN2k = + NZ^ , (31) 

C, = Zn\S„ (32) 

respectively. Note that the mapping /c 1— > A: is an involution of Z^r, i.e. k = k. Therefore 
it can be viewed as an action of the group Z2 on the space Z^- From this perspective. 
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Sg and Cg are nothing but the sets of points belonging to one-element and two-element 
orbits of the above action, respectively. The terms in the Hamiltonian involving pairs 
(fc, k) of fermions, in fact, depend only on the orbit. The XY Hamiltonian can be written 
accordingly as 



k&Sn 



H\^\g)=2j\ 
where 



COS 27r 



a + k 
1^ 



9 



kec„ 



(33) 



e N Ck 



e N c- 



(34) 



and h.y{k) is an hermitian operator on given by 

■ cos (27r^) -g i-f sin {2ir^) 



hy{k) 



-i7sin (27r^i±^) - cos (27r^) + ^ , 



(35) 



The factor 1/2 in front of the pair terms in fl55]) derives from the identity Cjhr^(k)Ck = 
Clh^{k)Ck, that expresses the fact that the various terms depend only on the orbit 
they belong to. 

Let us first focus on fermion pairs. For each k G Cg, h^{k) can be written as 



h^i{k) = —7 sin I 27r 



a + k 



N 



cos I 271 



a + k 



N 



9 



a 



(36) 



thus can be thought as a vector in the yz plane of the internal space of the pair, and 
is diagonalized (i.e. rotated up to the z direction) by a unitary rotation along x, 

RM)h^{k)RMy = ha\ (37) 

with /i G M and 



R 



x\Pk) 



■dk , 

exp ( -lyt^ 



cos 



-1 sm 



-1 sm Y 



cos 



(3^ 



cos Okcr^ + sinOka^ and R-JOkjcr^ R2 



x\yk 



it 



By recalling that Rx{9k)o-^Rx{ 
cosOkC^ — sinOkC^, and by requiring that the terms proportional to vanish, one 
obtains 

27r(a + k) 



. 2'K{a + k) 
7 sm — cos Ok 



cos ■ 



N 



9 



sin ^i- = 0. 



For each pair (/c, fc), there are two possible solutions that differ by vr. 



Ot — 6k 



STT, 



01 = -Ok + STT, 



SG{0,1}, 



where 



Ok = arctan 



7 sin {2n^) 
9 - cos (27r^) 



71 7T' 

"'-2' 2. 



(39) 



(40) 



(41) 
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Figure 1. Bogoliubov rotation for fermion pairs. 



and 



. 27r(a + fc) . 
- 7 sm — sm 9^^ + 

27r(a + k) 
cos — g 

N 



2n(a + k) 

cos ITr 9 

N 



cos 9t 



cosOl (l +tan^ 61) 



(42) 



The unitary transformation Rx{Ol) applied to Ck defines a new vector of fermion 
operators 

RM)Ck, keC„ SG{0,1}, (43) 



Bk,. 



where B^^ is related to B^^i by the relation 

Bk,i = i?x(7r)Sfc,o, k e C,. (44) 

See figure [H The fermion operators bk and 6^ are the Bogoliubov operators, and Rx is 
the Bogoliubov transformation for fermion pairs. By noting that 



cos^^ = (-1)^(1+ tan^^^) 



-1/2 



(45) 



for each pair of momenta one gets 

H!fl = Clh,{k)Ck = {-irel'\9) BlyBk,s, k e C„ se {0, 1}, (46) 
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where is the dispersion relation for fermion pairs 

27r(a + k) 



^k\9) = sgn 



cos ■ 



N 



cos ■ 



2n{a + k) 



. 2 27r(a + A;) 
+ ^2 51^2^— — ^. (47) 



Here sgn is the sign function, sgnx = x/\x\ for x ^ 0, and sgnO = 0. We stress that for 
each k E Cg the Bogohubov rotation is defined independently on the other pairs, and 
so the sign of the dispersion relation can be chosen in a completely arbitrary way pair 
by pair. It is not difficult to show that the unitary operator on the Fock space JF_(C^) 
corresponding to a Bogoliubov rotation Rx{0), 



bk = Uk{eiyckUk{ei) 

reads 

Uk{0) = exp 
Its action on the Hamiltonian is 



CfcCr + Cj^Ck, k G Cg 



:-iYe 



((>) 



k e Cg. 



(48) 
(49) 

(50) 



Observe that since Kk are quadratic with respect of creation and annihilation 
operators they commute with the parity operator fll4p . 

v = Ukie)vul{e), kec„ (51) 

and this means that the Bogoliubov transformation for fermion pairs preserves the parity 
sector. Finally, according to fl44p one gets the relation 

Uk{9l) = Vk-kUk{el), withVk-k = Uk{7r), keC,. (52) 

Note that the unitary operator V^^ can be decomposed in the form 

^kk = ^kkCkCk (53) 

where Ck and Sf^^ are respectively the charge conjugation and the swapping operator 



CkCkCf, — Cj^, 

^kk'^kSf^j^ = — iCfc, k G Cg, 

whose explicit expressions are 

Ck = exp (^i|(cfc + 4)) ' 
-^fcfc = exp (^i|(cicfc + c|cfe)) , keCi, 



(54a) 
(546) 



(55 a) 
(556) 

Consider now the case of single fermions, k G S^. The set depends both on the 
parity sector and on the parity of N. For even one gets 

'N' 

^ if ^ = -1, 

(56 a) 

if ^ = +1, 
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while, for odd, 

[{|01} if £>=-!, 

One can look at single fermions as a degenerate case of Bogoliubov pairs. Indeed, 
Equation fj^Tj) reduces to 

tanefc = 0, keSg, (57) 

whose solutions are given by Of. = svr, with s G {0, 1}. Therefore, in this case we are 
free to choose between two possible unitary transformation: the identity and the charge 
conjugation, 

Uk = {Ckf = l, or Uk = Ck, keS,. (58) 

Note that, if charge conjugation is chosen, parity is not preserved; rather the two parity 
sectors are swapped by the Bogoliubov transformation, 

V = -CkVCl (59) 

Finally, note that for single fermions the dispersion relation (H7|) reduces to 

ei'\g) = cos (stt^] - g, ke S„ (60) 



sin ( 27r^ 1=0, k e S,. (61) 



smce 

sin I ZTT — 

N 

In conclusion, the total Bogoliubov transformation that diagonalizes the 
Hamiltonian (!33l) has the form 

UB{g,r,Q,s)= H U,{ek){V,-,r l[{C,y^ , (62) 

where 

s = (sfc) e {0, 1}^, with Sk = Sfc, (63) 

and Cg/Z2 denotes that, in the case of coupled fermions, one must consider only one 
element for each pair (orbit of Z2). Due to the constraint in 0631) . the Bogoliubov 
unitary transformation has a gauge freedom represented by the arbitrary choice of a 
binary vector of length \Sg\ + \Cg\/2. 

Note that the anti-commutation relations are preserved by the Bogoliubov 
transformation, while the parity sectors are swapped according to 

V = {-ly^^^Usig, 7; Q, s) V Usig, 7; 0, sY, (64) 

where 

\s\e = \ss,\ = Sk. (65) 
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Therefore, one obtains the final expression of the diagonahzed Hamiltonian 



7 



2J 5^ (-1^5^(^)144 

k&N 



2 ' 



(66) 



where 



^=(-1)1^1'^^, (67) 

which depends on an arbitrary vector s G {0, Ijl'^el+I'^el/^^ that generates s by the relation 
Sk = sj, = Sk- Note that the physical part of H^^"^ acts on the sector of parity q. 



2.3. XY ground state: vacua competition 

One can use the gauge freedom of the Bogoliubov transformation ( l62l) in the following 
convenient way. Let s = s[g) be a function of the intensity of the magnetic field g, such 
that (— l)^fc(f) e\^\g) > for every g gM. From P7|) this means that 



sgn 



2TT(a + k) 
cos — g 



that is 



Sk[g) = - - - sgn 



N 



27r(a + A;) 

cos 77 g 

N 



(68) 



(69) 



Note that since Sk{g) = Sk{g), the above solution is consistent with the constraint (1631) 
of s. Therefore, the diagonahzed expression of the XY Hamiltonian reads 



2 ) ^^(9)' 



Qig) 



\\s{g)\, 



Q. (70) 



H\^\g) = 2J \e\:'\g)\ 

With this choice one has that in each parity sector the lowest energy state is the one 
with zero fermions {vacuum state) whose energy density / J is given by 



-^vac 



fceZiv 



g — cos 



2TTk 



1 2 



+ 7^ sin 



27rk 
1^ 



(71a) 



Eit) 



-E 



fceZiv 



g — cos 



2T^k 



TT 



+ 7^ sin 



^ + -1. (7U) 



Note, however, that a condition must be satisfied: the Bogoliubov vacuum state is a 
physical state, provided that it has the right parity Qig). Were this not the case, the 
projection Pg(g) would automatically rule it out. 

Let us look at the function Q[g) more closely. For even we have from (156 ap 
and (IMI) 



\^{g)\e = (s[0](^) + S[^](^)) 



1 + ^ sgn(l -g)-^ sgn(l + g) ) 5^ _i. 



(72) 
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whence 

Q{g) = {-l)\'\^Q = sgn(l - g^) + (73) 

that is 

Q{g) = sgn M - ^^9^) {N even) . (74) 



For odd we have from ( 15661) and ( 1691) 

|s(fi')le = S[o](5f)5g_i + sr]v-ii(5f)5g,+i 



11 1 



whence 



Q{g) = -sgn{l + Qg) (iVodd). (76) 

Since the vacuum state has holes, its parity is (—1)^, and it is a physical state 
only if 

q{9) = (-1)^- (77) 

Equation (iTTl) is satisfied for arbitrary q when g G (—1,1), while it is true only for 
Q = (—1)^ for g < —1, and g = +1 for g > 1. 

Therefore, for g G (—1,1), in the various regions of magnetic field g the ground 
state is alternatively given by one of the two vacua with energy f l71ap -( l7T6j) . We call 



this mechanism vacua competition between the two parity sectors. See Fig. [2l 

For < — 1 the vacuum state with ^ = — 1 for N even [g = +1 for odd) is not 
physical, because it has the wrong parity Q = Q = —{ — l)^ , and it is ruled out from the 
competition by the projection Pg. Analogously, for > 1 the vacuum state with g = —1 
for both A^ even and odd is ruled out. However, it is not difficult to prove that the 
energy of the unphysical vacuum when \g\ > 1 is always larger than the physical one. 
Therefore, as far as one is interested in the ground state, the ground state is the result 
of the vacua competition in the whole range (7 G M. Not so for the first excited level, 
which is the energy of the "losing" vacuum only in the range (—1, 1), while outside it is 
the lowest 1-fermion energy level above the losing vacuum. 

More generally, from (1741) and (176|) it easily follows that the whole spectrum is 
given for g & (—1, 1) by the union of the spectra of eigenstates with an even number of 
Bogoliubov fermions {g = (—1)^) of both Hamiltonians H^f^ with g = ±1. On the other 
hand, outside the above interval, the spectrum is given by the eigenstates of H^^^ {H^ 



7 



with an even (odd) number of Bogoliubov particles, where g = (—1)^ for g < —1 and 
g = +1 for g > 1. The intersection points between the vacua energy densities depend in 
general on the number of spins A^; however, independently of A^, the difference between 
the two energy densities. 



A^ 

meZ2jv 



V N J 
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A^ = 4 



E 




-0.6 
-0.8 
-1.0 
-1.2 



0.5 



IJO 



g 



-1.4 



N = b 



E 




9 



-1.4 



Figure 2. Vacua competition for = 4 and 5 spins: for both cases the dashed line 
corresponds to EiZc and the solid one to sltc ■ 



always vanishes at (7 = ±a/1 — 7^ (see figure [3]). Indeed one has 

1 



1 =F a/1 — 7^ cos ' 



/"Km\ 



and 



E (-1) (^) = ( i_e..(Ar+i)/^ ) 



0. 



(79) 



(80) 



From figure[3]on can also observe that for finite size systems the vacua intersection points 
present discontinuities of the first derivative, as will be explicitly shown in section HI In 
that section we will also focus on the points g = ±1 which are two interesting values of 
the magnetic field for this class of Hamiltonians, since they will be shown to represent 
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^diff 



vac 



N = A 



N = 5 




1.5 g 



Figure 3. Difference between energy densities (in unit J) of the two vacua at 7 = ^ 
for even numbers of spins, iV = 4, 5, 6. For all N, E^^^ vanishes at g = i/l — 7^. 



the finite-size forerunners of the quantum phase transition points (in the thermodynamic 
hmit). 

3. The XX model 



The XX model (7 = 0) is known as the isotropic model since the interaction between 
nearest neighbours spins along x and y axis is characterized by the same coefficient in 
the Hamiltonian ([T]): 



In this case equation 



reduces to 
a + k 

cos zvr 



9 



2 ' 



a 



(82) 



From this follows that the Fourier transformed XX Hamiltonian is already diagonal and 
the last term characterizing coupled fermions in Equation (l28l) vanishes for all k. In 
other words in the XX model we are only dealing with single fermions, Sg = Z^r, and 
the Bogoliubov transformation ( 162|) reduces to 

UB{g;Q,s)= H Cl\ (83) 



where now s G {0, 1}^ is an unconstrained binary string of length A^. This yields 
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Sk 



COS I 271 



a + k 



9 



Ckcl 



Pn 



(84) 



with Q = (— l)l*lf). In particular, if Sfc = the Bogohubov transformation associates 
JW fermions to Bogohubov fermions, while if = 1 it transforms JW fermions into 
Bogohubov antifermions, or holes. 



3.1. The energy spectrum 



As already emphasized at the end of section 12.21 the energy spectrum does not depend 
on the choice of the gauge s of the unitary Bogohubov transformation. If s = equation 
becomes 



(f) 



cos I 2n 



a + k 
1^ 



9 



c{ck 



15) 



The spectrum of the above Hamiltonian, and in particular its ground state energy has 
been studied in [22] . We quickly summarize the main results and show how they derive 
from vacua competition. The energy density/ J of the vacuum state does not depend on 
the parity g and on the size N 

a + k^ 



A9) = ^ E 



g — cos I 27T 



N 



9- 



(86) 



On the other hand, if we add one fermion of momentum k the energy density reads 



E['\k,g) = E^ 



g — cos I 2tt 



a + k 



a 



1 + Q 



(87) 



where g = —( — 1)^ is the parity of the 1-particle sector. In figure H] we represent the 
single particle energy spectra corresponding to = 8 and 9 sites (representative of an 
even/odd number of spins, respectively). The different lines are parametrized by /c G Z^r 
and one notes the presence of degeneracies in both cases. Since we are interested in the 
ground state of the system, we focus on the lowest energy levels and consider the values 
assumed by the function cos[27r {a + k)/N] in the four possible cases (iV even or odd and 
a = or 1/2 mod A^), as shown in Fig. [5l Notice that these results can be described in 
terms of regular polygons inscribed in a circle of unit radius, see figure [61 

From figure [5] one obtains the values of k that minimize the energy per site; in the 
1-particle sector one has 

A^ even =^ a = mod N ^ k = T^l , 

^ ' ^ r88) 
A^ odd ^ a = h mod N ^k - ^^'^ 



Similarly, in the 2-particle sector the energy is minimum for 

A^even a = 1 mod A^ ^ A;2} = { [f - l] , [f ] } , 

A^odd a = 0modN^{ki,k2}={[^],[^]}. 



(89) 
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N even {N = 8) 




N odd (A^ = 9) 




Figure 4. Single particle energy densities Ei{k,g) (solid lines) and vacuum energy 
densities £'vac(ff) (dashed lines). Different lines correspond to different k € 
according to ((57|) . 
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N = 9 Q = i mod N 




2-Kk I vr 
1^ + TV 



Figure 5. Plot of cos[27r(a + k)/N], k e Zjv for /V = 8,9 and a ee 0, 1/2 (modiV). 



It turns out that the general expression of the lowest energy levels in the different 
n-particle sectors does not depend on the parity of A^. For n fermions, one gets 

2 sin(n7r/iV) 



2n\ 



(90) 



N J N sin(7r/iV) ' 

In figure [7] we plot the lowest energy levels corresponding to < n < for N = 8 sites. 
The intersections of levels corresponding to n and n + 1 fermions (starting from n = 0) 
define the points of level crossing Qc, where an excited level and the ground state are 
interchanged. The analytic expression of the critical points is easily obtained by the 
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cos (a + A;)) 




2f (a + A;) 




^ cos 



(f(a + fc)) 



Figure 6. (Color online) Left: plot of cos[27r(a + k)/N], k e with a = (mod A^) 
(dashed line) and q = 1/2 (modA^) (solid line), for = 8. Right: geometrical 
description of cos[27r(a + k)/N] for = 8. 



condition Ef^'ig^) = E^l'Kg^). We find 



9c[n) 



-1 ' 



n+l 



1 + 2E 



m=l 



/7rm\ 

Kir) 



_ sm{mr/N) - sin[(n + 1)tt/N] 
~ sin(7r/Ar) ' 

As a consequence, the ground-state energy density is 

2 sin(n7r/A^) 



< n < iV- 1. 



(91) 
(92) 

ivs.n(./iv) ~r (^3) 

with < n < and where we stipulated that gc{—l) = — oo and gdN) = +oo. Thus, 
for g G {gc{n — l),gc{n)), the ground state contains n JW fermions. Note that gdO) = — 1 
and gc{N — 1) = +1, independently of A^. 

We will now derive the ground state energy density starting from the same choice 
of the Bogoliubov transform made for the XY Model fITU]) that particularizes to 



with g e (^gdn - l),gc{n) 



H^,'\g) = 2jJ2 



cos I 27r 



a + k 



clck 



2 ) ^^(9)' 



The ground state is then the winner of the vacua competition between 



-E 

k&N 

N ^ 



cos 



COS 



A^ 
27rk 



g 



(94) 
(95) 

(96 a) 
(966) 



(see Fig. [8]). The points of level crossing (19T]) are given by those values of the magnetic 
field that satisfy the following equation 



0. 



(97) 
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Figure 7. Lowest energy levels E"^"^{g) for different number of fermions n; the 
intersection between the energy levels corresponding to n and n + I fermions (starting 
from 71 = 0) are the points of level crossing (dots). 
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Figure 8. The ground state energy density of the XX model {N = 8) is given by the 
competition between the vacua energy densities E'vac (dashed line) and E^tc (solid 
line). 
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Figure 9. Geometrical representation of Equation ([57)1 when = 8; the thick red 
line is the magnetic field g. 



Consider the regular polygon inscribed in a circle of unit radius in figure [9l it is 
a geometrical representation of the function cos (^) for k G Zat. When \g\ > 1 one 
immediately gets E^^[g) = 0, whereas for \g\ < 1 the key idea is to consider the 
intervals on the x axis limited by the dashed vertical lines, represented in figure [U for 
each interval one can write the explicit expression for the vacua difference ( |97l) . For 
example when g G [cos (-^) , l] one gets 



N 
1 

N 



2(1-^7) 



g — cos 



nk 



Ilk 

q — cos I — 



(98) 



from which follows that in this interval E^^{g) = for = 1. Similarly when 



9 e 



cos 



(■m+l)7r 
N 



,COS^ 



for < m < — 1, one gets that E^!^{g) = when 
Tik 



cos 



k=l 



N 



-9c[m) 



<m< N -1, (99) 



where gdiTi) are the points of level crossing fl9Tl) [for n = one gets g = 1, in agreement 
with ( 1971) ]. By considering the symmetry gdn) = —g^N — 1 — n), one immediately 
sees that the level crossing points have the same analytic expression of the intersection 
points between the two vacua, for n = N — 1 — m. 
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4. Thermodynamic hmit and quantum phase transitions 

Quantum phase transitions in the XY model 

In this section we will show that in finite size systems one can find the forerunners of 
the points of quantum phase transition. These points are characterized by the presence 
of large values of the second derivative of the ground state energy density, that is then 
amplified and becomes a singularity in the thermodynamic limit. 

As observed in section [2^31 the first derivative of the ground state energy evaluated 
at the intersection points between the two vacua is not continuous and for finite size 
systems the second derivatives diverges at these points; however we will show that 
these singularities vanish when — oo. Consider for example the level crossing at 
g = a/1 — 7^; the difference between the first derivatives of the two vacuum energies is 
given by the derivative of fl78l) : 



v/Wcos(^ + ^) 
1 



(100) 



1- v/Wcos(2^). 

When the number of spins is odd, for each k E ^Ln there is a given k = k + ^ such 
that cos (^) = — cos + (see Fig. [TU]), and the last equation becomes 



dE'dltt / . N 1 ^ 



diff , ^ 1 ^,2 



dg 'J N ^ 



r 



2v/W+4^"g'' cos(^)v/W 



y l-cos(^)(l-y^ 



;(101) 



from the symmetries of the function cos (^) one gets that the last expression is 
strictly greater than zero. From this it follows that the second derivative of the vacua 
energy difference diverges for all finite N at g = y^l — 7^, and the same argument 
can be extended to all intersection points between the two vacua. The case A^ even, 
see figure El is analogous, as one can see by noting that the polygon corresponding 
to a = i(modA^) is rotated by an angle (or in other words, it associates to each 
momentum k ^ k = k + ^^j^)- 

Summarizing, for finite size systems the second derivative of the energy density of 
the ground state diverges at the intersection points of the two vacua; on the other hand 
in the thermodynamic limit this divergence is suppressed. Indeed, in the limit N ^ 00 
equation fllOOp becomes: 



where f{x) is given by: 



m = - — , , (103) 

1 — a/1 — 7 cos(xj 
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Figure 10. Geometrical description for cos { ^'^^'^'^^ ^ for iV = 9, similar to figure [6j 
the two polygons belong to the parity sectors a = (modA^) (dashed line) and 
a = 1/2 (modiV) (solid line) 



Expanding in Taylor series f{x 
dE: 



n/N) one gets 



?diS 
^vac 



r 



=/'(0) = 0, 



N 



oo. 



^ . V ^ , ; , . v^y - ™- (104) 

dg \' 'J n^/T^ 
This means that the singularities of the second derivative of the ground state vanish 
in the thermodynamic limit; in other words, the forerunners of the quantum phase 
transition are not related to finite-size level crossings of the ground state. In this section 
we will show that they are related to the level crossings between the unphysical vacuum 
and the losing physical vacuum where single Bogoliubov fermions sit. 

Consider the explicit expressions of the vacua energies corresponding to the four 
possible cases given by the parity of and the two parity sectors 

(i) iVeven, f? = -l,5,= {[0], [f]}, 



(ii) N odd, g 



m 




;(105a) 



(1056) 



(105c) 



fceZiv 
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(iv) iVodd, f? = +l,5,= {[^]}, 

vac jY )/ 

(105d) 

Observe that the absolute values in the previous expressions correspond to the cosines 
evaluated at single fermion momenta Sg\ at these values of the magnetic field the first 
derivative of energy is not continuous (see figure [2]) and the second derivative has terms 
proportional to the Dirac delta functions 5{g±l). However, remember that the vacuum 
in case (i) becomes unphysical as soon as \g\ > 1, so that aX g = ±1 there is a level 
crossings between physical and unphysical states. The same phenomenon happens to 
the vacuum in case (ii) a.t g = 1, and to the vacuum in case (iv) at g = —1. On the 
other hand, one can observe that for finite size chains, for both even and odd N, the 
ground state is smooth at = ±1, in other words the ground state, which coincides with 
the winning vacuum state, does not have any singularities at these points. However, it 
can be shown that the second derivative of the ground state energy at g = ±1 scales 
as — log . 

Consider for example the case of an even number of spins A^. In this case the 
ground state belongs to the parity sector with g = +1, without singularities. Figure [Tl] 
displays d'^Eitc/dg'^ for A^ = 6,24,54; at 5? = ±1 it scales like — logA^. Indeed when 



cos 



2TTk 



IT 

N 



7^ sin^ 



27rk 



TX 
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Figure 12. Second derivative of the vacuum energy density at 17 = 1 in the parity 
sector Q = +1 for N even, from 18 to 320 (dashed Hne): it scales as — logA^ (solid 
line). 



(? = 1, by deriving ( 17161) one has 



^ 4; Isin + #)| [1 + 7^ + cos + f ) (y - 1 



(l + cos(^ + f)) 



3/2 



3/2 



1 

77r 



3 + log 



7 



2^3/2 



:i+7 

-1 and N even {g 



— ^logiV, 
77r 



(106) 
1) and g = ±1 



for — > 00, as shown in figure [121 The cases : 
and odd {g = =i=l) are analogous. 

The quantum phase transition is forerun by the losing vacuum whose second 
derivative contains a Dirac delta function, at the transition between physical and 
unphysical states. When tends to infinity, as we will now show, the difference between 
the two vacua at g = ±1 tends to zero and quantum phase transition forerunners 
approach the ground state, building up singularities at logarithmic rates. Indeed, at 
= ±1 from Equation (!78|) one has: 



2'Kk 



TX 

N 



(107) 



where 



/±(a;) = Y (±1 — cosx) + 7^ sin x. (108) 

In the thermodynamic limit, by applying the same technique used in (I102p . equation 
(11071) becomes 

'.27r 



e: 



diS, 



1 

±1 ~ 

^ 2n 



fix) 



dx 
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Figure 13. Difference between the two vacua energy densities at g = ±1: exact result 
(dotted line) and asymptotic approximation of order (solid line). 
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27r 2!iV2 



2tv 



f"{x)dx 



IT 



77, 



oo, 



(109) 



2N^ 

where we used the equahty /4(0) = —7. See figure [131 In figure [H] we display the low 
energy part of the spectrum (thin lines) and the energy density of the two vacua (thick 
lines): at (7 = ±1 the ground state is the winning vacuum that has no singularities, the 
first excited level coincide with the losing vacuum for (7 G (— 1, 1). Its second derivative 
diverges at g = ±1, forerunning the quantum phase transitions. Observe that they 
are at the transition between a physical state, which coincides with the first excited 
level, and an unphysical state, which does not corresponds to any physical level: for 
l^fl > 1 the losing vacuum is unphysical. Summarizing, we identify as forerunners of 
the quantum phase transition those points of the losing vacuum energy density whose 
second derivative diverges. These points are associated to single Bogoliubov fermions 
and belong to the crossing between the first excited level and the unphysical vacuum 
for finite size systems. When iV — *• 00 they approach the ground state as N~^. 



4-2. Quantum phase transitions in the XX model 

As observed in Section [3] the XX model (7 = 0) is characterized by the only presence 
of single fermions, and the absence of Bogoliubov pairs. As a result, all points 
g = cos ^ ^^(^+^) j (in both parity sectors) with k G Tjn can be considered quantum 
phase transitions forerunners. See fl96ap -f l96fej) and compare with fll05ap -f ri05(il) . Indeed, 
the second derivative of the vacua energy density contains a Dirac delta function at 
these points and, apart from g = ±1, they all belong to the first excited level like in 
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Figure 14. The thin Unes represent the low energy part of the spectrum of the XY 
chain, and the thick solid and dashed lines refer to the two vacua energy densities, Mac 
and iSvac , respectively: these vacua alternatively coincide with the ground state and 
the first excited state for \g\ < 1. When g = ±1 one vacuum energy is the ground state 
energy, while the other one docs not corresponds to any physical level. The transition 
points are the forerunners of the quantum phase transition. 



the XY model, see Fig. [15] (we will focus on g = ±1 at the end of this section). In the 
thermodynamic limit these points forerunning the quantum phase transition approach 
the ground state, becoming critical points. Consider for example gi = cos (^) 7^ ±1; 
the energy difference between the vacua is now given by 



1 

N 
1 

N 



COS 



COS 



27li 



27li 



COS 



COS 



27rk 



COS 



27li 



— COS 



27rk 



7T 

N 



2txI 



IT 

N 



(110) 



XY model on the circle 

N = A E 



26 




g 



-1.5 -1.0 -0.5 




9 



Figure 15. The thin lines represent the low energy part of the spectrum of the XX 
chain with iV = 4 and N = 5 spins; the solid and dashed tick lines refer the two vacua 
energy densities in the parity sectors with g ^ —1 and g = +1, respectively. The 
forerunners of the (continuous) quantum phase transition points arc indicated with 
bold points; they are given by gk = cos {^ff-) , k G Zjv- 



By using the same technique of the previous section one gets 



-^vac i9i 



where 
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Figure 16. Difference between vacuum energies at £ = 3 versus N (daslied line), and 
its asymptotic approximation (solid line). 

From the symmetries of f{x) and its derivatives, it follows that equation fillip becomes 

See figure [161 Therefore, in the thermodynamic limit the forerunners of the quantum 
phase transition in the isotropic XX model approach the ground state faster than the 
ones of the XY model (with 7 7^ 0). Compare figures fT3] and fT6l 

As shown in figure [T5l the intersection points of the two vacua (which coincide with 
the level crossing points gi{n) discussed in section [3]) are characterized by a discontinuity 
of the first derivative for finite size chains. By deriving the energy difference flOTjl . one 
can show that the discontinuity of the first derivative at the points of level crossing scales 
like 1/A^; therefore in the thermodynamic limit the divergence of the second derivative 
vanishes, as for the XY Hamiltonian with 7 7^ 0. 

Let us finally consider the points g = ±1: on one hand they are level crossing points 
{g = ±^1 — 7^, 7 = 0), on the other hand, following the same criterion introduced for 
the XY model, they can be considered as forerunners of quantum phase transitions: 
what happens in this particular case is that these points belong to the ground state 
already for finite A^. Another crucial difference between the anisotropic case and the 
XX model is that, since all Bogoliubov fermions are single, there are A^ + 1 points 
forerunning the quantum phase transition. Thus in the A^ ^ 00 limit they densely 
fill the interval [—1, 1] of g and yield, as one expects [1], a continuous quantum phase 
transition in this interval. 
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Conclusions 

In this paper we analyzed the XY model with periodic boundary conditions. Being 
interested in finite size systems, we did not neglect the boundary term which derives 
from the Jordan- Wigner transformation. In order to diagonalize the Hamiltonian we 
deformed the discrete Fourier transform with a local gauge coefficient depending on the 
parity of spins down, anti-parallel to the magnetic field. We then showed that in the 
Fourier space there are two classes of fermions, single and coupled ones; this distinction 
is crucial in order to determine the Bogoliubov transformation, which is also gauge 
dependent. From the expression of the diagonalized Hamiltonian we reinterpreted the 
ground state and the first excited level of the system as given by a competition between 
the vacuum energies of the two parity sectors. We finally introduced a criterion to 
find those values of the magnetic field that can be considered forerunning quantum 
phase transitions in the thermodynamic limit. They are associated to single Bogoliubov 
fermions and to the level crossings between physical and unphysical states. 

There is considerable interest in the study of entanglement for quantum spin chains, 
both in view of applications and because of their fundamental interest. See, for example, 
the results concerning the XX chain [T6l 1251122]. Future activity will focus on the study 
of the properties of the multipartite entanglement of the ground state in terms of the 
distribution of bipartite entanglement [26l EI] and on the investigation of the possible 
connections with quantum phase transitions in the thermodynamic limit. 
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